And quiet flows the supersolid 4 He 
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A superfluid having atomic scale superflow of a hexagonal lattice of vortex and antivortex fila- 
ments, described by a single macroscopic wave function is presented as a supersolid. As superfluid 
4 He is pressurized, at a first order transition, rotons (atomic scale current circulation, a vortex loop) 
not only condense but also expand and fuse into hexagonal or other complex superflow patterns. The 
vortex core contains an excess density of non-condensate atoms. Further, a Kelvin (m = 0, necklace) 
mode condenses in the vortex filaments. It results in a 3D atom density wave of hep symmetry. In 
our theory, superfluid phase stiffness, rather than atom localization, imitates a solid like rigidity. 
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Supersolid [1,2] is a quantum crystal that exhibits non 
classical moment of inertia [3], a type of superfluid re- 
sponse. Decades of efforts [4], following early theoretical 
suggestions, have culminated in a recent striking observa- 
tion of non classical moment of inertia, by Kim and Chan 
[5]. This is yet another jewel in the crown of condensed 
4 He , an elusively simple one component bose system; the 
deeper one searches and digs, the more surprises are in 
store. This work has excited a renewed interest [6,7] in 
the quantum many body theory of supersolid. It is likely 
to open new directions and find new phenomena in ex- 
periments; it also offers a good opportunity to bring the 
field of cold atom BEC, boson Mott insulators etc., closer 
to supersolid 4 He . A correct theoretical understanding 
of the supersolid phenomena will throw new light into 
old experimental quantum anomalies in the superfluid- 
solid helium interface [8], such as crystallization waves, 
Kapitza resistance etc. 

The supersolid predicted by Andreev, Lifshitz [1] and 
Chester [2] has a reference crystal of localized atoms; 
large amplitude quantum fluctuations create ground 
state vacancies or defects which undrgo condensation. 
In our mechanism superfluid is the reference system and 
solid like rigidity emerges from the superfluid stiffness in 
a fundamental way. A spatially periodic, atom scale su- 
perfluid flow develops spontaneously in the ground state 
leading to a solid like response, in addition to superflu- 
idity. We call it a roton fusion(ROFU) mechanism, as 
the atom scale flow pattern arises spontaneously from a 
condensation and a complex fusion of real rotons (vortex 
loop of atomic dimensions), as we pressurize superfluid 
4 He . ROFU mechanism is one collective way of releas- 
ing the kinetic energy frustration, that is present in a 
solid with localized atoms. It is a kind of spontaneous 
generation of staggered Abrikosov vortex lattice, created 
by hydrostatic pressure. Hydrostatic pressure, a scalar, 
mimics an internal staggered magnetic field. 

Roton minimum of superfluid He has been viewed 
[9], for quite some time as a soft mode that drives a 
superfiuid-solid phase transition, as rotons have wave 
vectors close to the reciprocal lattice vectors of the hep 
4 He crystal. A new and crucial ingredient in our theory is 



to use the non-trivial inner structure of roton, an atomic 
size quantized vortex loop. In our theory rotons do more 
than condensation; they expand and fuse into ordered 
vortex and antivortex filaments (atom scale thickness) 
and loose their identity. 
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FIG. 1. Figure 1. a) Region of roton minima in the 
(k x , k y , 0) plane. Qi, Q2 and Q3 are Fourier component vec- 
tors of the macroscopic condensate wave function (^(r)). b) 
Unit cell of the hexagonal vortex-antivortex array (tubes); it 
has a 2D atom density wave of hexagonal symmetry c) Kelvin 
(azimuthal quantum number m = 0, necklace) mode with a 
finite ±Q Z condenses in the vortex in a staggered fashion. It 
results in an atom density wave of hep lattice. Lattice sites 
(not position of localized atoms) are denoted by shaded circle. 

In superfluid 4 He , nearly 90 % of the atoms are out 
side the zero momentum condensate. These large frac- 
tion of atoms with finite momenta appear as spatially 
uniform quantum fluctuation of the superfluid vacuum. 
They are part of the vacuum and the ground state is 
a 100 % superfluid. When vortices and antivortices are 
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spontaneously created in the ground state, a part of the 
non-condensate fraction gets piled up in the vortex core 
region, creating a net atom density wave of hep sym- 
metry. The self consistent potential that maintains this 
redistribution arises from the underlying supcrfluid stiff- 
ness. In this sense supcrfluid stiffness gives a solid like 
rigidity to the 3D atom density wave. A fragile 2% su- 
pcrfluid condensate(as measured by Kim and Chan), by 
a complex maneuver of its own flow, orchestrates the 3D 
atom density wave and a rich 'lattice dynamics'. Periodic 
superflow leads to a periodic modulation of quantum fluc- 
tuation and non condensate component. 

In what follows we develop some heuristic pictures, fol- 
lowed by a Bogoliubov theory and then discuss some con- 
sequences. 

Kim and Chan have discovered non-classical moment 
of inertia in pressurized solid 4 He ; they observe strong 
thermal hysteresis, small variation of supersolid transi- 
tion temperature T ss for a pressures range, 25 to 60 bars, 
and an anomalous increase of T ss with addition of 3 He 
impurity. More importantly, the authors indicate that 
the superfluid response may not be due to zero point va- 
cancies or defects or interfaces. That is, it is likely one 
is dealing with a solid with integer number of atoms per 
unit cell and yet exhibiting a supcrfluid property. 

While building our theory we remember that at the 
supcrfluid-solid first order transition at T = 0, inter 
4 He atom distance decreases only by a small, ~ 3%, 
(a density decrease ~ 10%). This means that locally the 
quantum solid is no more crowded than the quantum liq- 
uid is. Thus any generalized rigidity that emerges on 
both sides is likely to have the same quantum charac- 
ter, namely some kind of phase stiffness (ODLRO) aris- 
ing from local coherent number fluctuations. Further, 
the superfluid solid coexistence line in the P-T plane 
has nearly zero slope for T < 0.2 K indicating that the 
two phases have the same entropy per mole, according to 
Clausius-Clapeyron equation ^f- = ff-Epf- This means 
that density of states of low energy bosonic quasi parti- 
cles is nearly the same for superfluid and supersolid at 
the coexistence line, indicating a possible deep connec- 
tion between the stiffness of a supcrfluid and rigidity of 
a supersolid. 

Now we present a Bogoliubov theory to illustrate our 
ROFU mechanism. Bogoliubov theory is a mean field 
theory that works well for weakly interacting bosons. 
However, for a proper choice two body pseudo poten- 
tial, it captures qualitative features and some quanti- 
tative features, including roton spectrum. The model 
Hamiltonian is: 

H = ]T( efe - n)bl b k + ^ ]T VilHK-A^'+ci (!) 

Here €k = ^ is the kinetic energy of 4 He atoms with 
mass M, chemical potential \i, created by operators &t's. 



V(q) is the two body, effective or pseudo potential and 
ft is the volume of the system. 

In the uniform superfluid phase, Bogoliubov theory 
starts by replacing the zero momentum operator 6 — > 

nfie^Ni, a classical expectation value, with no as a 
condensate fraction. The Bogoliubov quasi particle spec- 
trum of the uniform supcrfluid state has the well known 
form 

«"» - \/(w) 2 + (w) 2 "» v w < 2 » 

where no is the zero momentum condensate fraction. 
When the potential V(k) is negative for a range of k, the 
spectrum has a roton minimum; an example is Bruckner 
pseudo potential, V(k) « V a ^%^-. with a w 2.1 A and 
V «15if. 

Landau suggested roton as the first excitation that has 
associated with it a rotational velocity flow (and hence 
the name roton). Feynman viewed it as a single atom 
motion, but dressed by a back flow, the net effect being 
a quantized vortex loop of atomic dimensions. Recent 
variational study and wave packet analysis [10] confirms 
this and shows roton wave packet as a ball of current 
disturbance, and little density variation. 

Within the Bogoliubov theory, the roton minimum (fig- 
ure la) reaches zero energy (become completely soft), 
when V(k) — — ^r^TTST - This indicates an instability 
of the uniform superfluid solution. We expect a roton 
condensation and reorganization of the ground state. A 
complex current flow associated with a roton, rather than 
a simple density variation, made us wonder what will 
be the consequence of a roton condensation. It became 
clear that a primary consequence of roton condensation 
will be formation of a current density wave, not atom 
density wave. This resulted in our roton fusion hypoth- 
esis. According to this hypothesis the condensed rotons 
expand and fuse into a hexagonal array of vortices and 
antivortices. Our hypothesis is naturally influenced by 
the known hep structure of solid 4 He . 

To test our hypothesis we have compared the energy of 
uniform superfluid state with an off diagonal long range 
order (ODLRO) ansatz that contains a hexagonal array 
of vortices and antivortices. After some struggle we found 
a simple and elegant ansatz - in phase superposition of 
three plane waves generates the desired vortex structure: 

(*(r)) = ^o(r) ~ 4e i *(e iQl " + e lQ2 ' r + e lQ3 ' r ) (3) 

It also gives a lower energy than the uniform supcr- 
fluid state, for a choice of V(k). Here Qi = ^f(;4p0), 

Q.2,3 = 2s(^|,±l), Qi + Q 2 + Q 3 = 0, 'a' is the lat- 
tice parameter of the hexagonal lattice. Here, no and 
'a' are variational parameters. The above ansatz leads 
to variation of atom density only in the xy-plane, a 2D 
atom density wave. This results in a supersolid behaviour 
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along the x and y directions and superfluid in z-direction 
- it is an anisotropic supersolid, a type of quantum liquid 
crystal. Later we will improve it to get a 3D atom den- 
sity wave of hep symmetry. We can find more complex 
flow pattern (braided and knotted vortices) leading to an 
atom density wave of hep symmetry using more Q's etc. 

One of the triangular sub lattices of the hexagonal lat- 
tice contains vortices and the other antivortices (figure 
lb). It is easy to show that equation (3) has asymptotic 
form 

(<]>(r)) <~ x + iy (vortex), x — iy (antivortex) (4) 

respectively around the vortex and antivortex filaments 
of the two sub lattices. 

In addition to global U(l) symmetry (overall phase of 
(^(r))) and translational symmetry there is a discrete 
symmetry breaking arising from P&T violation. That is, 
another degenerate solution, not connected by a global 
phase rotation of equation (5) is obtained by the replace- 
ment Qi — ► Qi in equation (3); cquivalently, by an 
interchange of vortices and antivortices. 

According to equation (4) the condensate fraction van- 
ishes quadratically as we approach the core of line vor- 
tices. Further, the fraction of the total particles con- 
densed in our macroscopic wave function is < 1. Where 
are the rest of non-condensed particles ? In the Bogoli- 
ubov theory of uniform superfluid, these are the parti- 
cles that carry non zero momenta and they appear as a 
spatially uniform quantum fluctuations of the zero mo- 
mentum condensate. They influence, for example the su- 
perfluid stiffness, through the spectrum; otherwise these 
finite momentum components are not visible in the low 
energy dynamics of the irrotational superfluid. 

In our ROFU, a spatial density variation of the non- 
condensate fraction is induced. In particular they appear 
as extra atom density in the normal core region of our 
ordered vortices, where the condensate density vanishes. 
This is somewhat counter intuitive, as one expects a de- 
pletion of fluid density at the vortex core. It is known 
from an early work of Fetter [11] that within Bogoliubov 
theory, in the vortex core region the local atom density 
increases compared to the uniform background atom den- 
sity. This is a remarkable non-local quantum effect, ap- 
pearing within the Bogoliubov theory, as explained by 
Fetter. 

We elaborate this point further, by comparing a refer- 
ence Bose Einstein condensed (BEC) state 

*bec ~ n( eiQir + e * Q2 ' r + e * Q3 ' r ) ( 5 ) 

i 

with the N-particle projected wave function of Bogoli- 
ubov theory for our ROFU solution: 

*N[r 4 ] = ^x(r P i,r P2 )x(rp3,rp4) \{ y pn-i, r PN ) 

p 

(6) 



Summation over permutation P symmetrises the 
wave function. The pair function x(ri,r2) = 
77(ri,r 2 )V'o [ ri 2 1 " 2 ] • When the pair function r?(ri,r 2 ) = 
1, equation (6) becomes the same as BEC (equation 5). 
In the BEC state, atom density, by construction, van- 
ishes as we approach the core of the line vortices. In the 
Bogoliubov wave function, equation (6), the pair func- 
tion rj makes an important difference and creates a a pile 
up of atom density through pushing some of the non- 
condensate fraction to core regions. The pair function 77 
represents, a repulsive correlation induced between any 
two particles in the Bogoliubov theory. It shows how 
pairs of particle are pushed in and out of the condensate 
by the Bogoliubov process (&Q 1 )(6Q 2 )^fc t _ k+Ql+Q2 etc., 
in a space dependent fashion. 

We want to make certain remarks about the nature 
of ODLRO in our supersolid phase. Our primary order 
is a long range order in momentum space, an ODLRO. 
Bosons condense in a single macroscopic wave function 
V'o(r) given by equation (5). It has three Fourier coeffi- 
cients which are not independent: (6q 1 )= (&q 2 ) = (^q 3 ) 

= nje 1 ^: they have same amplitude and phase. 

The spatial variation of atom density implied by our 
ODLRO should not be thought of as a diagonal long 
range order (DLRO). A spatial periodic ordering of a 
small fraction of the non-condensate fraction is induced 
self consistently in our theory. For example in our Bo- 
goliubov factorization we get an anomalous term such 
as ^Qi^Q2^k^k-Qi+Q 2 ' -^ density wave of wave vec- 
tor (Q x — Q 2 ) induced by the above term has its origin 
in single particle condensation. In principle we could 
generate such an anomalous term through non- vanishing 

averages such as, Sfc(^k^k+Qi-Q 2 )- This wm ^ c an m " 
dependent diagonal long range order (DLRO) parameter. 

The solution we have discussed so far, equation (3), 
has no density variation along z-axis. In real 4 He we 
expect maximum amount of atom density to be concen- 
trated in the normal core region. Thus it has its natural 
tendency for local spatial order, arising from short range 
repulsions. We view this as a Kelvin (azimuthal quan- 
tum number m = 0, necklace or sausage ) mode condenses 
(figure lc) at wave vectors Q z = in the vortex fila- 
ments. We modify equation (3) and generate a periodic 
modulation of the vortex core size, along the vortex fila- 
ment in a staggered fashion, 

(*(r)) ~ (x + iy)(l + ecosQ z z) and 

~ {x -iy)(l-e cos Q z z) (7) 

close to core of the vortices of the two sub lattices. In 
view of the short range interaction of the model Hamil- 
tonian, the total energy is reduced further for small e, 
another variational parameter. This vortex core modula- 
tion leads to an atom density wave of hep symmetry and 
a 3D solid like rigidity. 
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We can choose our model parameters of Hamiltonian 
(equation 1) such that our hexagonal lattice contains an 
average of 2 atoms per unit cell or close to it. 

From our solution it follows that the superfluid density 
is mostly concentrated outside the vortex core, which is 
the interstitial region of the hep atom density wave (fig- 
ure lb and lc). 

So far we have sketched a Bogoliubov theory, a ROFU 
solution and some key features of the ground state. Our 
theory is far from rigorous and complete. However, the 
physical argument for a spontaneous generation of mi- 
croscopic circulating ground state current is compelling. 
Bogoliubov (mean field) theory, in view of its non pertur- 
bative character is capable of finding possible new phases 
in dense liquid 4 He . 

In what follows we show how our theory qualitatively 
explains salient features of Kim-Chan's results. Further, 
we briefly sketch interesting consequences that also fol- 
low; some of them are very unique to our theory. 

In our ROFU solution we discussed a 2D solid and a 3D 
solid. It is likely that as a function of pressure there is a 
small region in the pressure-temperature diagram where 
the 2D solid intervenes, as shown in figure 2. 



vortices and thereby increase the supersolid-solid transi- 
tion temperature. Both are consistent with Kim-Chan's 
observations. The anomalously low critical velocity, for 
destruction of the non classical moment of inertia (super 
solidity), observed by Kim and Chan is likely to arise 
from a collective depinning of the ground state vortices 
from their weak self consistent potential. 

A key prediction of our theory is the presence of atom 
scale circulation. One important consequence of this is 
on lattice dynamics. The P & T violation in the ground 
state and the first order dynamics of vortices makes the 
lattice dynamics different from the classical hep solid. We 
get splitting of degeneracies etc.; some of the anomalous 
modes can be viewed as coupled Kelvin modes. We will 
discuss them in a future publication. 

In summary, our supersolid is a superfluid in disguise 
with ineresting consequences. 

It is a pleasure to thank N. Kumar for an early (1971) 
and insightful introduction to Leggett's paper [3] and a 
recent discussion; Dung-Hai Lee for alerting me of Kim 
and Chan's work [5] ; Moses Chan for a brief discussion at 
Hong Kong, implying 'no room for vacancy' superfluidity 
in their results. 
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FIG. 2. Schematic phase diagram for 4 Hc in the pressure 
temperature plane. An intermediate region of anisotropic 
supersolid (supersolid in x, y direction and superfluid in 
z-direction) is suggested. 

We can qualitatively discuss some features of finite 
temperature phase transitions that are unique to ROFU 
solution. Thermally produced vortices in the ROFU con- 
densate will proliferate, depin and melt the underlying 
vortex lattice. Disappearance of ODLRO will result in a 
solid of localized atoms. As the process is a melting of 
a vortex lattice, i) phase transition will be a first order 
one and ii) small traces of 3 He atoms could help pin the 
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